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1. Introduction 

A few observational and/or experimental results have dramatically pushed forward the 
research program on gravity as those from the radio-metric Doppler tracking received from 
the Pioneer 10 and 11 spacecrafts when the space vehicles were at heliocentric distances 
between 20 and 70 Astronomical Units (AU). These data have conclusively demonstrated the 
presence of an anomalous, tiny and blue-shifted frequency drift that changes smoothly at a 
rate of~6xl0^^ Hz s^^. Those signals, if interpreted as a gravitational pull of the Sun on 
each Pioneer vehicle, translates into a deceleration of ap = (8.74 ± 1.33) x 10"^" m s~^. This 
sxmward acceleration appears to be a violation of Newton's inverse-square law of gravitation, 
and is referred to as the Pioneer anomaly, the nature of which remains still elusive to unveil. 
Within the theoretical framework of nonlinear electrodynamics (NLED) in what follows we 
wiU address this astrodynamical puzzle, which over the last fifteen years has challenged 
in a fundamental basis our understanding of gravitational physics. To this goal we will 
first, and briefly, review the history of the Pioneers 10 and 11 missions. Then a synopsis 
of currently available Lagrangian formulations of NLED is given. And finally, we present 
our solution of this enigma by invoking a special class of NLED theories featuring a proper 
description of electromagnetic phenomena taking place in environments where the strength 
of the (electro)magnetic fields in the backgroimd is decidedly low. 
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2 What is the problem: The Pioneer anomaly 

In this short voyage to the Pioneer 10 and 11 missions our main guide will be the 
comprehensive and richly documented recent review on the Pioneer Anomaly by [Turyshev, 
S. G. & Toth, V. T. (2010). Living Rev. Rel. 13 (2010) 4. arXiv:1001.3686, v2, gr-qc] from 
which we retake some ideas and references. (The attentive readers are kinldy addressed to 
this invaluable article). 

The Pioneer 10 and 11 spacecrafts were the first two man-made space vehicles designed to 
explore the outer solar system. The trajectories of the spaceships were projected to passage 
nearby Jupiter during 1972-1973 having as objectives to conduct exploratory investigation 
of the interplanetary medium beyond the orbit of Mars, the nature of the asteroid belt, the 
environmental and atmospheric characteristics of Jupiter and Saturn (for Pioneer 11), and to 
investigate the solar system beyond the orbit of the Jovian planetFl 




The Pioneer missions were the first space probes to adventure over the asteroid belt, heading 
for close-up observations of the gaseous giant planets, and for performing in situ studies of the 
physical properties of the interplanetary medium in the outer solar system. The design of their 
missions was guided by the simplicity, having a powerful rocket-launching system to push the 
spacecrafts on an hyperbolic trajectory aimed directly at Jupiter, which the spacecrafts were 
expected to fly-by approximately 21 months after laimch (see Fig. 1). 

By the late 1960's, the aerospace engineering technology available to the designers of the 
Pioneer missions made it no longer practical to use solar panels for operating a spacecraft 
at large distances, as for instance that of Jupiter. A cause of this, a built-in nuclear power 
system, in the form of radioisotope thermoelectric generators (RTGs) powered by ^^Fu, was 
chosen as the means to provide electrical power to the spaceship. As even this was relatively 
new technology at the time the missions were designed, the power subsystem was suitably 
over-engineered, being the unique design requirement to have a completely functional space 
probe capable of performing all planned scientific tasks by running only three (out of four) 



The entire design of these spacecrafts and their science missions was characterized by such 
conservative engineering, and for sure it was responsible for both the exceptional longevity 
of the two spacecrafts and their ability to deliver science results which by far exceeded the 
expectations of their designers. 

The original plan envisioned a primary mission of two to three years in duration. 
Nevertheless, following its encounter with Jupiter, Pioneer 10 remained functional for over 
30 years. Meanwhile, Pioneer 11, though not as long lived as its engineering-copy craft, 
successfully navigated a path across the solar system for another encounter with Saturn, 
offering the first close-up observations of the ringed planet. After the encounters with Jupiter 
and Saturn (for Pioneer 11, see Fig. 1), the space ships followed, near the plane of the ecliptic, 
hyperbolic orbits of escape heading to opposite sides of the solar system, continuing their 
extended missions. The spacecrafts explored the outer regions of the solar system, studying 
energetic particles from the Sun (solar wind), and cosmic rays entering our neighborhood 
in the Milky Way. (Their cousin spacecrafts, the Voyagerl and 2, that where launched 
contemporarily, studied in the beginning of their mission, the interplanetary space, what 
resulted in a very accurate mapping of the interplanetary magnetic field and its strength, as 
one can see in Fig. 2 below). 

^See details on the Pioneer missions at http : / /w ww . nasa . gov/ cent ers/ames /miss ions /| 
[archive /pioneer .html Be awared that another member of Pioneer spacecrafts family. Pioneer 6, 
remained operational for more than 35 years after launch. 




RTGs. 
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Fig. 1. Ecliptic pole view of the spacecrafts Pioneer 10 ar\d 11 interplanetary trajectories (see 
also the trajectories of the vehicles Voyager 1 and 2). Credit: 

http : / /www . nasa . gov/ cent ers/ames /miss ions /archive /pioneer . html 



In virtue of a combination of many factors, the Pioneers were excellent space sondes for 
pursuing experiments of high precision celestial mechanics. This includes the presence of 
a coherent mode transceiver on board, the attitude control (spin-stabilized, with a minimum 
number of attitude correction maneuvers using thrusters), the design of the power system (the 
RTGs being on extended booms aided the stability of the craft and also helped in reducing 
thermal effects), and Doppler tracking of high precision (with the accuracy of post-fit Doppler 
residuals at the level of mHz). The exceptional built-in sensitivity to acceleration of the 
Pioneer 10 and 11 spacecrafts naturally allowed them to reach a level of accuracy of ~ 10^^" 
m/s^. The result was one of the most precise spacecraft navigations in deep space since the 
early days of space exploration. That is the great legacy of the Pioneer missions. 
After having had a brief accoimting of the Pioneers missions, one can proceed to review our 
current understanding of nonlinear electrodynamics and to settle down the foundations for its 
use in the search for a solution to the Pioneer anomaly. In this Section we shall briefly review 
the theoretical foundations of some theories of NLED, focusing essentially on the fundamental 
prediction concerning the way photons propagate through a vacuum space permeated by 
electromagnetic (EM) fields: The fact that photons travel along the effective metric, and not 
over the geometry in the background. It is this peculiar feature what makes the photon to 
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"feel" itself being acted upon by a force, and consequently to imdergo acceleration. |^ In our 
understanding, such effect is responsible for the drift in frequency undergone by the photon. 
Next we will show that any NLED, independently of the specific form of its Lagrangian, brings 
in such a frequency shift. And in our view, it is such acceleration what can account for the 
Pioneer anomaly. 

3. Some Lagrangian formulations of nonlinear electrodynamics 

To start with, it is worth to recall that according to quantum electrodynamics (QED: see 
[Delphenich j2003 2006} for a complete review on NLED and QED) a vacuum has nonlinear 



properties (Heisenberg & Euler 1936; Schwinger 1951) which affect the photon propagation. 
A noticeable advance in the realization of this theoretical prediction has been provided by 
[Burke, Field, Horton-Smith , etal., 1997), who demonstrated experimentally that the inelastic 
scattering of laser photons by gamma-rays in a background magnetic ield is definitely a 
nonlinear phenomenon. The propagation of photons in NLED has been examined by several 
authors [Bialynicka-Birula & Bialynicki-Birula, 1970; Garcia & Plebanski, 1989; Dittrich & Gies, 
1998; De Lorenci, Klippert, Novello, etal, 2000; Denisov, Denisova & Svertilov, 2001a, 2001b, 
Denisov & Svertilov, 2003]. In the geometric optics approximation, it was shown by [Novello, 
De Lorenci, Salim & etal., 2000; Novello & Salim, 2001], that when the photon propagation 
is identified with the propagation of discontinuities of the EM field in a nonlinear regime, a 
remarkable feature appears: The discontinuities propagate along null geodesies of an effective 
geometry which depends on the EM field on the background. This means that the NLED 
interaction can be geometrized. An immediate consequence of this NLED property is the 
prediction of the phenomenon dubbed as photon acceleration, which is nothing else than a 
shift in the frequency of any photon traveling over background electromagnetic fields. The 
consequences of this formalism are examined next. 

3.1 Heisenberg-Euler approach 

The Heisenberg-Euler Lagrangian for nonlinear electrodynamics (up to order 2 in the 
trimcated infinite series of terms involving F) has the form Heisenberg & Euler jl936| 



LH-E = -^F + af' + ^G2, (1) 

where F = F,,^Ff^, with Ff,v = d^^Ay - dyA^„ and G = {nafi-ysF^f^F^^ = -4E ■ B, with greek 
index running (0, 1, 2, 3), while a and ^ are arbitrary constants. 

When this Lagrangian is used to describe the photon dynamics the equations for the EM 
field in vacuum coincide in their form with the equations for a continuum medium in which 
the electric permittivity and magnetic permeability tensors and y.^^ are fimctions of the 
electric and magnetic fields determined by some observer represented by its 4-vector velocity 
y [Denisov, Denisova & Svertilov, 2001a, 2001b; Denisov & Svertilov, 2003; Mosquera Cuesta 
& Salim, 2004a, 2004b]. The attentive reader must notice that this first order approximation 
is valid only for B-fields smaller than Bq = = 4.41 x 10^^ G (Schwinger 's critical B-field 
|Schwinger|^1951^ ). In curved spacetime, these equations are written as 

^Because of the special theory of relativity constraints regarding the propagation of any perturbation, 
it becomes clear that such effect must manifest itself as a change in one or both of their physical properties: 
its frequency or its wavelength. Hence, through the Pioneer spacecrafts radio Doppler tracking we might 
be observing the effect on the photon frequency. 
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Yl. 

• c 



(2) 
(3) 
(4) 



Here, the vertical bars subscript "||" stands for covariant derivative and ^"^f"^ is the 
antisymmetric Levi-Civita tensor. 

The 4-vectors representing the electric and magnetic fields are defined as usual in terms of the 
electric and magnetic fields tensor F^v and polarization tensor Pj,y 



c 



Hu = P, 



t'" c ' 



c ' c ■ 

where the dual tensor X*^ is defined as X*^, = jtjfivafiX"'^ , for any antisymmetric second-order 
tensor X^^. 

The meaning of the vectors D'' and H'' comes from the Lagrangian of the EM field, and in the 
vacuum case they are given by 



(5) 



(6) 



where the permeability and tensors are given as 



1 + 
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2a 
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2a 
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7a 



7a 



(7) 



(8) 



(9) 



In these expressions a is the EM coupling constant {<^ = Yc~ TS?)- ^^'^ tensor h^y is the metric 
induced in the reference frame perpendicular to the observers determined by the vector field 

Meanwhile, as we are assioming that = 0, then one gets 



P P 457rB2 P 
and /i„fl = f(?!aR. The scalars e and \i can be read directly from Eqs.|8|l9| as 



la 



457rB2 



B2 



(10) 



(11) 



Applying conditions l |62) and > 63 (derived in the Appendix) to the field equations when £* = 
0, we obtain the constraints efe^mk^' = and bf'kfi = and the following equations for the 
discontinuity fields and ba : 
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ya 



(12) 



(13) 



Isolating the discontinuity field from | (T2| |, substituting in equation (13) , and expressing the 
products of the completely anti-symmetric tensors flv^y^V'^''^^' ^ terms of delta fimctions 
|Stephani|(2004^ , we obtain 



+ 



k^ + 



fiB^b^jk^kPf 



(^j^lj.hn^B'^k^vf'^V^ =Q. (14) 



This expression is already squared in k^, but still has an unknown fe^ term. To get rid of it, one 
multiplies by B^, to take advantage of the EM wave polarization dependence. By noting that if 
B* feff =0 one obtains the dispersion relation by separating out the A:''fc^' term, what remains is the 
(-) effective metric. Similarly, if Ba-fo" ^ 0, one simply divides by B^fcT so that by factoring out 
, what results is the (+) effective metric. For the case Bfffe* = 0, one obtains the standard 
dispersion relation 



whereas for the case B^ fe" ^ 0, the result is 



0. 



(15) 



1 + 



^B2 



xf, _¥_ 



}l }l- ^B2 j " }l 
where (') stands for and we have defined 

7ci 



+ 



457rB2' 



and 



+ 



kaka = Q, 



Bt' 
I B-TB^, 1 1/2 



(16) 



(17) 



as the unit 4-vector along the B-field direction. 

From the above expressions we can read the effective metric g"^ and g'^, where the labels 
"+" and "-" refers to extraordinary and ordinary polarized rays, respectively. Then, we need 
the covariant form of the metric tensor, which is obtained from the expression defining the 
inverse metric gfivg""^ — ^fi- So that one gets from one side 



S}tv — Spy 



and from the other 



(18) 



+ 



u'B 



i^)(l + ^ 



S}IV 
Vf,Vy 



. A''B 
1 + ^ + 



^B2 



fi-pB^ 



1 + 



-1 



gpy 



Uv ■ 



(19) 
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The function ^ can be expressed in terms of the magnetic permeability of the vacuum, and 
is given as 

^=2fl-iV (20) 
Thus equation (19) indicates that the photon propagates on an effective metric. 
3.2 Born-lnfeld theory 

The propagation of light can also be viewed within the framework of the Born-lnfeld 
Lagrangian. Such theory is inspired in the special theory of relativity, and indeed it 
incorporates the principle of relativity in its construction, since the fact that nothing can travel 
faster than light in a vacuum is used as a guide to establishing the existence of an upper limit 
for the strength of electric fields around an isolated charge, an electron for instance. Such 
charge is then forced to have a characteristic size [Born & Infeld l |1934| . The Lagrangian then 
reads 



(21) 



As in this particular case, the Lagrangian is a fimctional of the invariant F, i.e., L = L(F), but 
not of the invariant G = Bp,Ef', the study of the NLED effects turns out to be simpler (here 

again we suppose E = 0). In the equation above, b = = = = 9.8 x 10^^ e.s.u. 

" SI 

In order to derive the effective metric that can be deduced from the B-I Lagrangian, one has 
therefore to work out, as suggested in the Appendix, the derivatives of the Lagrangian with 
respect to the invariant F. The first and second derivatives then reads 

^F = — ^77172 a'^d Lff=-— 372- (22) 

The L(F) B-I Lagrangian produces, according to Eq.l|70| in the Appendix, an effective 
contravariant metric given as 

C = , ~],U2 S'' + , "'3/2 [^"^ + '"'1 • (23) 
4(1 + 1^) 2fo2(l+F) 

Both the tensor hfn, and the vector in this equation were defined earlier (see Eqs.lj9j and | [I6| 
above). 

Because the geodesic equation of the discontinuity (that defines the effective metric, see the 
Appendix) is conformal invariant, one can multiply this last equation by the conformal factor 

/ p\3/2 

-4(1 + 52) to obtain 

S^f={^ + ^)s''-^-^[h'"+l''n- (24) 

Then, by noting that 

F = Ff,yF^"■' = -2{E^ - B^) , (25) 
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and recalling our assumption £ = 0, then one obtains F = 2B^. Therefore, the effective metric 
reads 



JiV 



2B^ \ 2B^ 



g'^-^w'+^fn, (26) 



5eff - 1,2 JS 1,2 

or equivalently 

C = r + ^-^v"v^-^-^pr. ill) 

As one can check, this effective metric is a functional of the background metric g'"', the 
4- vector velocity field of the inertial observers V" , and the spatial configuration (orientation 
V) and strength of the B-field. 

Thus the covariant form of the backgroimd metric can be obtained by computing the inverse 
of the effective metric g-^g just derived. With the definition of the inverse metric g^^gfl = <5'a, 
the covariant form of the effective metric then reads 



- Sliv (2B2/fo2 + 1) + (2B2/fo2 + 1) ' 



vi'hich is the result that we were looking for. The terms additional to the background metric 
gfiv characterize any effective metric. 

3.3 Pagels-Tomboulis Abelian theory 

In 1978, the Pagels-Tomboulis nonlinear Lagrangian for electrodynamics appeared as an 
effective model of an Abelian theory introduced to describe a perturbative gluodynamics 
model. It was intended to investigate the non trivial aspects of quantum-chromodynamics 
(QCD ) like the asymptotic freedom and quark confinement [Pagels & Tomboulis| jl978|) . In 
fact, Pagels and Tomboulis argued that: 

"since in asymptotically free Yang-Mills theories the quantum ground state is not controlled by 
perturbation theory, there is no a priori reason to believe that individual orbits corresponding to minima 
of the classical action dominate the Euclidean functional integral. " 

In view of this drawback, of the at the time understanding of groimd states in quantum 
theory, they decided to examine and classify the vacua of the quantum gauge theory. To 
this goal, they introduced an effective action in which the gauge field coupling constant g is 

r pi pi7 }n' -1 

replaced by the effective coupling g{t) • T = In I -^^^^^j — I ■ The vacua of this model correspond 
to paramagnetism and perfect paramagnetism, for which the gauge field is F^^j, = 0, and 
ferromagnetism, for which F^^yF" = A^, which implies the occurrence of spontaneous 

magnetization of the vacuum. They also found no evidence for instanton solutions to the 
quantum effective action. They solved the equations for a point classical source of color spin, 
which indicates that in the limit of spontaneous magnetization the infrared energy of the field 
becomes linearly divergent. This leads to bag formation, and to an electric Meissner effect 
confining the bag contents. 

This effective model for the low energy (3+1) QCD reduces, in the Abelian sector, to a 
nonlinear theory of electrodynamics whose density Lagrangian L{X,Y) is a functional of the 
invariants X = FjivF^^ and their dual Y = (FfiyF^^)*, having their equations of motion given 
by 

^This is the imprint that such theory describes nonlinear electrodynamics. 
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\/f,i-LxF^'" -LY*Ft"')=0, (29) 

where Lx=dL/dX and Ly = dL/dY. This equation is supplemented by the Faraday equation, 
i. e., the electromagnetic field tensor cyclic identity (which remains unchanged) 



V;,F,a + V,Fa,,+VaF,,^ 

In the case of a simple dependence on X, the equations of motion turn out to be KimzeH2008^ 
(here we put C = and 47 = — (A^)(*^^)/^ in the original Lagrangian given in Pagels & 
rTomboulis]p78) ) 



0. 



(30) 



as; 



(*-l)/2 



X, 



(31) 



where (5 is an dimensionless parameter and [A] = {anenergy scale). The value S = 1 yields the 
standard Maxwell electrodynamics. 

The energy-momentum tensor for this Lagrangian L(X) can be computed by following the 
standard recipe, which then gives 



Tuv — 



47r 



(32) 



while its trace turns out to be 



r = 



AS j 



X. 



(33) 



It can be shown 



Kxmze 



2008 I that the positivity of the = p component implies that i5 > 1/2. 



The Lagrangian \31) has been studied by [Kunze {2008^ for explaining the amplification of 



the primordial magnetic field in the Universe, being the analysis focus ed on three different 
regimes: 1) > E^, 2) ~ 0{E^), 3) <C B^. It has also been used by |Mosquera Cuesta & 
[Lambiase, (2009) to discuss both the origin of the baryon asymmetry in the tmiverse and the 
origin of primordial magnetic fields. More recently it has also been discussed in the review on 
" Primordial magneto-genesis" by Kandus l |2010| . 

Because the equation of motion (29) above, exhibits similar mathematical aspect as eq. (35) 
(reproduced in the Section), it appears clear that the Pagels and Tomboulis Lagrangian (31) 
leads also to an effective metric identical to that one given in equation (40), below. 



3.4 Novello-Perez Bergliaffa-Salim NLED 

More recently, Novello et al. 1 2004 1 (NPS) revisited the several general properties of nonlinear 
electrodynamics by assuming that the action for the electromagnetic field is that of Maxwell 
with an extra term, namel}|^ 



(34) 



where F = Fj^yF^" . 



^Notice that this Lagrangian is gauge invariant, and that hence charge conservation is guaranteed in 
this theory. 
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Physical motivations for bringing in this theory have been provided in Novello et al. (2004 1. 
Besides of those arguments, an equally unavoidable motivation comes from the introduction 
in the 1920's of both the Heisenberg-Euler and Born-Infeld nonlinear electrodynamics 
discussed above, which are valid in the regime of extremely high magnetic field strengths, i.e. 
near the Sch winger's limit. Both theories have been extensively investigated in the literature 
(see for instance Mosquera Cuesta & Salim j2004a[b| ; |Mosquera Cuesta et al. l |2006 and the 
long list of references therein). Since in nature non only such very strong magnetic fields exist, 
then it appears to be promising to investigate also those super weak field frontiers. From the 
conceptual point of view, this phenomenological action has the advantage that it involves only 
the electromagnetic field, and does not invoke entities that have not been observed (like scalar 
fields) and/ or speculative ideas (like higher-dimensions and brane worlds). 
At first, one notices that for high values of the field F, the dynamics resembles Maxwell's 
one except for small corrections associate to the parameter 7, while at low strengths of F it is 
the 1/F term that dominates. (Clearly, this term should dramatically affect, for instance, the 
photon-B field interaction in tntergalactic space, which is relevant to understand the solution 
to the Pioneer anomaly using NLED.). The consistency of this theory with observations, 
including the recovery of the well-stablished Coulomb law, was shown in [Novello et al.|(2004^ 
using the cosmic microwave radiation bound, and also after discussing the anomaly in the 
dynamics of Pioneer 10 spacecraft Mbelek et al. |2007|. Bo th analysis provide small enough 
values for the coupling constant 7 , Mosquera Cuesta pOlO l 



3.4.1 Photon dynamics in NPS NLED: Effective geometry 

Next we investigate the effects of nonlinearities in the evolution of EM waves in the vacuum 
permeated by backgroimd B-fields. An EM wave is described onwards as the surface of 
discontinuity of the EM field. Extremizing the Lagrangian L(F), with F(A^,), with respect 
to the potentials yields the following field equation Plebanski 1 1970^ 



Vy{LFF^'")=0, (35) 
where Vv defines the co variant derivative. Besides this, we have the EM field cyclic identity 



Taking t he discontir i uities of the field Eg. (35) one gets (all the definitions introduced here are 
given in 



Hadamard 



1903| )p] 



Lpf/k^+2LFFF''Pf,^Ft'^k,,=0, 
which together with the discontinuity of the Bianchi identity yields 



(37) 



/ff/j^T + fiii^fi + fh^oi = . (38) 
A scalar relation can be obtained if we contract this equation with k'^F'^^, which yields 



^Following Hadamard's method 



Hadamard 



1903|, the surface of discontinuity of the EM field is 



denoted by Z. The field is continuous when crossing E, while its first derivative presents a finite 



discontinuity. These properties are specified as follows: [F^ 



= 0, 



[f,,i,|a 



where the 



symbol [F^n'] j; = lim^'^o+ (/Is+iS ^ /Is-i) represents the discontinuity of the arbitrary function / through 
the surface Z. The tensor /j,„ is called the discontinuity of the field, k,\ = is the propagation vector, 
and the symbols "1" and "m" stand for partial and covariant derivatives. 
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(f*%^r''"+2F''V)fcF^v = 0. (39) 

It is straightforward to see that here we find two distinct solutions: a) when f^/a^ = 0, case 
in which such mode propagates along standard null geodesies, and b) when F"^ f^^ = %■ In 
the case a) it is important to notice that in the absence of charge currents, this discontinuity 
describe the prop agation of the wave f ront as determined by the field equation j35) , above. 
Thence, following Lichnerowicz. 1 1962 1 the quantity can be decomposed in terms of the 
propagation vector fca and a space-like vector (orthogonal to fca) that describes the wave 
polarization. Thus, only the light-ray having polarization and direction of propagation such 
that P^f'knap = will follow geodesies in Any other light-ray will propagate on the 
effective metric jio). Meanwhile, in this last case, we obtain from equations p7 and (|39l 



the propagation equation for the field discontinuities being given by Novello et al 



2000) 



kuh = . 



(40) 



effective metric 

This equation proves that photons propagate following a geodesic that is not that one on 
the background space-time, gf", but rather they follow the effective metric given by Eq.|40), 
which depends on the background field f", i. e., on the B-field. 

4. Understanding the Pioneer anomaly within NLED 
4.1 Astrodynamics of Pioneer 10 and 11: Input facts 

As pointed out above, since 1998 the JPL group have continuously reported an anomalous 
frequency shift derived from about ten years study of radio-metric data from Pioneer 
10: 03/01/1987-22/0 7/1998 [Anderson et al.| p98) . Pioneer 11: 05/01/1987-01/10/1990 
[Anderson e t 31.^1995). The group has also found a similar feature in the data from of Ulysses 
and Galileo spacecrafts Anderson et al.'fl998, 2002). The observed effect mimics a constant 



sunward acceleration acting on the spacecraft with magnitude 



flp = (8.74 ±1.33) X 10" 
and a steady frequency (v) drift 



cm s 



which equates to a "clock acceleration": 



: 6 X 10"^ Hz/s 



(41) 



(42) 



(43) 



dAv 

dt c 

where c represents tha speed of light in a vacuum, and t is the one way signal travel time. An 
independent analysis of the radio-metric Doppler tracking data from the Pioneer 10 spacecraft 
for the period 1987 - 1994 confirms the previous observations ^Markwardt p002j . In addition, 
by removing the spin-rate change contribution yields an apparent anomalous acceleration 
flp = (7.84 ± 0.01) X 10^^ cm s^^, of the same amount for iDoth Pioneer 10/11 Anderson et 
[aL||2002) ; [Abramovici & Vager jl986) . Besides, it has been noted that the magnitude of flp 
compares nicely to cHq, where Hg is the Hubble parameter today. 
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As stressed above, unlike other spacecrafts like the Voyagers and Cassini which are 
three-axis stabilized (hence, not well-suited for a precise reconstitution of trajectory because 
of numerous attitude controls), the Pioneer 10/11, Ulysses and the by-now destroyed Galileo 
are attitude-stabilized by spinning about an axis (parallel to the axis of the high-gain 
antenna) which permits precise acceleration estimations to the level of 10^^ cm s^^ (single 
measurement accuracy averaged over 5 days). Besides, because of the proximity of Ulysses 
and Galileo to the Sun, the data from both spacecrafts were strongly correlated to the solar 
radiation pressure unlike the data from the remote Pioneer 10/11. Let us point out that the 



motions of the four spacecrafts are modelled by general relativistic equations (see Anderson 
[et al.^ |2002l, section IV) including the perturbations from heavenly bodies as small as the 
large main-belt asteroids (the Sun, the Moon and the nine planets are treated as point masses). 
Nonetheless, the observed frequency shift remains unexplained Turyshev & Toth| l |2010) . 
Thenceforth, several proposals for dedicated missions to test the Pioneer anomaly are now 
under consideration [Pioneer Collaboration (2(305 1, in virtue of the dramatic implications of 
the Pioneer puzzle for the understanding of gravity. 



4.2 What has been done by other researchers 

In search for a possible origin of the anomalous blueshift, a niimber of gravitational and 
non-gravitational potential causes have been ruled out by Anderson et al. j2002) . According 
to the authors, none of these effects may explain ap and some are 3 orders of magnitude or 
more too small. The addition of a Yukawa force to the Newtonian law does not work ease. 
An additional acceleration is predicted by taking into account the Solar quadrupole moment 



Mbelek & Michalski| l |2002) . Although this entails a blueshift, it decreases like the inverse of 



the power four of the heliocentric radius, being of the order of ap only below 2.1 AU 
Meanwhile, the claim that the Modified Newtonian Dynamics (MOND) may explain ap in the 
strongly Newtonian limit of MOND Quevedo| f2005| ; IMilgron] (2001| |200 2) is not obvious at 
all. First, the fits to the rotational curves of spiral galaxies yield for the MOND acceleration 
constant flg ^ value eight times smaller than cHq Quevedo |2005^ . Second, the gravitational 
pulling of the Sun up to 100 AU still yields an acceleration greater than aq by at least three 
orders of magnitude, equating ag only at about 3000 AU. Hence, Newtonian dynamics up 
to general relativity corrections should apply to the spacecrafts. Otherwise, one would be 
inclined to conclude that MOND is ruled out by a laboratory experiment Milgron|p001 2002, 1. 
Now, any true Doppler shift would involve an accompanying acceleration, which would be in 
conflict with both the motions of planets and long-period comets Anderson et al. l |1995 ; lorio 
1 2006a b e I. 

Heretofore what we have learnt is that based on Einstein-Maxwell equations, the only other 
photon frequency shift that can be misinterpreted, at the solar system scale, with the Doppler 
shift is the gravitational frequency shift. In the weak field and low velocity limit, this would 
involve a time dependent gravitational potential instead of a spatial dependent one. Such 
proposals invoking the dark energy as the source of the time dependent gravitational potential 
have been suggested |Iorio| |2006a b c l; Tangen (2006). However, quintessence, like other 
fundamental scalar fields, has not yet been observed. 

In simmiary, prosaic explanations, non-gravitational forces and modified dynamics or new 
interaction (long or short range) force terms do not work'Mbele k & MichaIski|(2002| ; [Quevedo 



^2005 l; Milgron ( 2001 : 2002 |; Ranada ( 2003; 2005) . Gra vitational origin of the anomaly is rouled 



out by the precision of the planetary ephemeris (see Anderson et al. jl998^ , [lono| j2006a b c >, 
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and others Tangen ( 2006)) and the known boiinds on dark matter within the orbital radius of 
Uranus or Neptune Raflada]j2003[|2005^ ; [Miitmire & Matese|j2003f . 



4.3 What we are proposing to tackle the Pioneer anomaly 

By gathering together all the arguments reviewed above, one is led to the conclusion that the 
Pioneer anomaly does not seem to be related to the gravitational interaction [Anderson et aL] 
{1998} ; [lorio| |2006alblc| ; [Tangen[ p006) . If this is the case, what other of the currently known 



interactions in nature could afford a consistent imderstanding for the radio-metric Doppler 
tracking data from Pioneer spacecrafts? 

The right answer could be related to the fact that there are only two long range interactions 
known today: Gravity and electromagnetism. Therefore, what remains is the EM sector|^ 
Meanwhile, the possibility of an interaction of the EM signal with the solar wind leading to a 
change of the frequency of the EM signal is now rouled out (see |Anderson et aL] j2002| ). 
Indeed, it appears to be unescapable to conclude that what we are observing (measuring 
through the receivers) could be related to the equation of motion of the photon. In other 
words, the mounting evidence seems to converge to what could be happening to the photon 
during its propagation through the interplanetary space from the Pioneer 10/11 antennas to 
the receivers on Earth. 

It is timely, then, to recall that classical (Maxwell theory) or quantized (QED) linear 
electrodynamics does not allow for a change of the frequency of a photon during its 
propagation in a linear medium without invoking diffusion due to the interaction with the 
surroimding matter (hence a smear out of the image of the source). Moreover, for such a 
phenomenon to occur, one needs to consider a general and non trivial Lagrangian density 
L = L(F) for which its second derivative w.r.t. F: d'^L/dF^ = Lff ^ 0. Therefore, the Pioneer 
anomaly, if not an artifact, may be a result of NLED as we show below. Indeed, relation l )43| 
above translates, in covariant notation, into 

^V,fc/' = ^fc'', (44) 
dl 

where / is some affine parameter along a ray defined by = ^ (see'Fujii & Sasaki' /2006^). 
The latter equation departs from the classical electrodynamics one (see ^Landau & Litchiftz] 
jl970f , section 87) 

^Vvfc''=0 (45) 

and suggests the occurrence of the NLED effect dubbed photon acceleration. 
The concept of photon acceleration, which follows from the description of photon propagation 
in NLED, was discussed by [NoveUo & Salim |2001) , see also the book by ^Mendonga et aL] 



j2006| . Next we explain why the anomaly shows up in some situations and not others. (For 



experimental tests of NLED and further theoretical predictions see 'Mosquera Cuesta et al 
(20061; B urke et al.[p997) ; [LurLdstrom et aL[ j2006) ; [Lundin et al . (20061; Mar klimd & Shukla 
(2006f ). 

Therefore, the alternative that the Pioneer anomaly is not consequence of an actual change 
in the spacecraft velocity (see [Anderson et aL] j2002| . Section X) deserves to be investigated. 
Indeed, a direct interpretation of the observational data from the spacecrafts implies merely 
an anomalous time-dependent blueshift of the photons of the communication signals. On the 



'^Non-metric fields can also be regarded as gravitational fields and there is a lot of space for speculation. 
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other hand, in using a time dependent potential lorio 12006a b e); Tangen (2006) to explain the 
Pioneer 10/11 data one may be pointing out to the need of an effective metric for the photons. 
In fact, what is needed is just a time variation of the 4-momentum of the photon along its path. 
Thus the atomic energy levels would not be affected. Rather, only the motion of the photon 
being concerned. 

4.4 NLED at all distance scales: From cosmology down to astrodynamics in the Solar 
System 

Upon the collection of arguments presented above, it appears that all these requirements are 
achieved by considering that NLED is based on a Lagrangian density L(F) w hich includes 
terms depending nonlinearly on the invariant F = Fj,t, , withF = 2(B^c^ — £^) 



Novello et al. 



j2000f ; |Novello & Saliinlj2001f ; pebanski|fl970| , instead of the usual Lagrangian density L 



|F of classical electromagnetism in a vacuum. As stated above, we shall explore the effects 
of nonlinearities in the evolution of EM waves, which are envisioned onwards as the surface of 
discontinuity of the EM field. Therefore, as shown above, by extremizin g the Lagrangian w ith 
respect to the potentials A^, one obtains the EM field equation of motion Plebanski < 1970^ 



V4LfF'"')=0, (46) 
in which Vy represents the covariant derivative, and Lf = dL/dF. 

Recalling the discussion above, the dynamics of the photon propagation follows the equation 

r''"-4^Ff'«F/)fc,,fcv = 0. (47) 

which exhibits the fundamental feature of NLED, i.e., the effective metric. 
Then, by taking the derivative of the last expression, one arrives to 

r\/yk,=A(^Fl'l'F^'kuh) . (48) 



Eq.(j48| shows that the nonlinear Lagrangian introduces a term acting as a force accelerating 
the photon. This acceleration of any photon which is traversing over weak background 
electromagnetic fields in a vacuum is the new physical element that we argue hereafter would 
be responsible for the Pioneer anomaly. 

4.4.1 NLED photon acceleration: What Earth receivers are reading off radio signals from 
interplanetary spacecraft transponders - The case Pioneer anomaly 

If NLED is to play a significant role at the macroscopic scale, this should occur at the 
intermediary scales of clusters of galaxies or the interclusters medium, wherein most 
observations show that the magnetic fields are almost uniform (and of the same order of 
magnitud^, unlike the dipolar magnetic fields of the Sim and planets. However, galaxies 
are gravitationally bound systems, whereas the cosmic expansion is acting at the cluster of 
galaxies scale. Thus, the magnetic field (B) in clusters of galaxies (IGMF) depends on the 
cosmic time (B = Bga^^). So, the B that is relevant to this study is that of the local cluster 



^Next we show that the "acceleration" of photons predicted by NLED may account for the anomalous 
blueshift indicated by the Pioneer 10/11, Ulysses and Galileo spacecrafts. This will manifest itself as a new 
frequency shift for the EM waves, in addition to the Doppler shift (special relativity) and the gravitational 
an d cosmological redshi ft (general relativity), when both of them apply. 



*Fuju& Sasaki 2006 
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of gal axies Beck (2000). (As regard to the contribution of the CMB radiation see Riess etalj 
j2004| )p] Recently, \Meej ^2002 1 has speculated that the 2 fiG magnetic field he has observed 
within the local supercluster of galaxies in cells of sizes of about 100 kpc may extend all the 
way down to the Sun. We explore further this idea in the framework of NLED and show that 
it is capable to provide an explanation of the Pioneer anomaly from first principles. 
Relation l40l can be cast in the form 



(49) 



where V = F^^ky and h^ = hl' h^, . 

As E = 0, one can write, after averaging over the angular-dependence Bialynicka-Birula & 
|Bialynicki-Birula|(T970} : 



-l\m 



with ||fc|| = CO I c = 1nv I c. By inserting this relation in 49' yields 



Taking the (x*) derivative of Eq.|5T) we obtain 

2g^vk^[r)u + kn'[g^x, 



c2 ^ Lf 



c2 U 



(50) 



(51) 



(52) 



The cosmological expansion will be represented by g^y = fl2(j/)§-|,'"'^'''' , with a the scale factor, 
^{local) ^j^^ IqquI metric. So, Eq.' 52 1 yields: 

,2 



i] the conformal time, and g^y 

2g,,yk^k')\0 + 2 



g^,vkV = 



-FF 



Lf 



where the dot stands for partial derivative w.r.t. rj. Using Eqs.i51 

O^^LpF _ 1 fLO^pl^FF 



and 1 53 we obtairp' 



fc,<(fc'')|o = 



(53) 



(54) 



Lf 2 V c2 Lf y |o 

Now, f = — 4 (^) F, by recalling that (x a^^. Moreover, from the method of the effective 
metric, it can be shown that fc" does not vary with time in the first order approximation unlike 
||fc||[^Hence 



'The interclusters magnetic field is in any case by far small (10^' G) to add a measurable correction 
even to the cosmological redshift. As for the contribution of the cosmic microwave background (CMB), 
not only it is too weak but also, the CMB is pure radiation (F = 0), whereas we are interested in the case 
of a background magnetic field with no significant electric field counter-part, i.e., £ = 0. 

^"By removing the NLED extra term from Eq.|49|, this reduces it to = so that the photons 

would just see the local background metric. 

^^Given a background metric gjan as a result of NLED effects photons follow geodesic paths with 

respect to the effective metric (or any one conformal to it) gfj^'^ — gfw — 4^F^f F^i, (see 



Novello et al. 



(2000|[Novello & Salim| {2001} ). Thus, following the usual analysis on the gravitational frequency shift 
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CO I CO 

c \ c 



By inserting relation |55} in |54| , and then expanding and arranging, one finds 

'a\ Q + 2FQf 



Q 



(55) 



(56) 



where we have set Q = F ^ and Qf = dQ/dF. 

At present cosmological time (t), and for a duration very short as compared to the universe 
age, Eq.(|56| reduces to 



■H, 



Q + 2FQf 
1-Q 



(57) 



where (v is the photon frequency t-derivative). v 7^ if and only if a) the NLED contribution 
is non-null, i.e., Lpp 7^ 0, and b) F depends on time. 

4.4.2 NLED alia NPS as explanation of the Pioneer puzzle 

The explicit form of this general nonlinear Lagrangian (which simulates the effect of dark 
energy in |Novello et al.| f2004 1) reads 



1^ 7 
I'+F' 



4 F« ' 



(58) 



where n is a strictly positive integer. From Eqs.i 56[58) , the time variation of the photon 
frequency, due to interaction with very weak B{t) fields, reads 



4n7„ - (2fi + l)F"+i 



(59) 



(F"+i + 4ny„) [F«+i + in{n + 2)^„] ' 

with A„ = 4HoH(n + 1). Notice that 7,, should be negative in order to guarantee that the 
Lagrangian is boimd from below (see ^Landau &: Lifchiftz] ^1970) , sections 27 and 93), 7,; = 
— (B„c)^("+^\ Also, it is worth noticing that Eq.i59i in the nearly-zero field limit (B — > 0) 
would reduce to 



n + 1 
n + 2 



(60) 



which implies a blueshift. 



4.5 Discussion and conclusion 

We stress that the NLED is a imiversal theory for the electromagnetic field, with 7„=i = 7 in 
Eq.ljSSf being a universal constant. The value of 7 was fixed in Novello et al. j2004| by using 



but with gjiif-'^' replacing gjiy, one gets k°c = u}q/\J gg^l^^ (see _ 

where Vg = cvq/Itt denotes the photon frequency in flat Minkowski spacenme. 

cosmological redshift (subsequent to the time dependence of the curvature), the variation of /c" with time 
can be neglected in the first order approximation, since F'^''F^^ = F^Fj^q = in the case of a zero electric 
field. 



Landau & Lifchiftz 



19701, section 88), 
rtius, discarding the 
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the CMB constraintp^ Indeed, in the standard model of cosmology the CMB is well described 
by Maxwell theory which is likely to give a good accoimt of the magnetic fields in galaxies too. 
Notwithstanding, the processes at the origin of both the seed magnetic field in the interclusters 
medium and clusters of galaxies are not yet clearly understood. Hence, if NLED is to play a 
significant role at the macroscopic scale, this should occur at the intermediary scales of clusters 
of galaxies or the interclusters medium. 

Then, considering the possibility that the NLED correction terms described above come into 
play at these scales, one gets the following ordering of strengths: Byniverse ^ Bintercluster ^ 
< Bgaiaxy Turning back to the Pioneer anomaly, a good accordance is obtained 



Bi < B, 



Cluster 



with Blsc = 1.7 7 Bi (case n 
0.008 ± 0.002 ffGlmessetal 



2004 



l)|Riess etal. 1 2004 >. Thus, setting Bi = 1 ItI^^*, one finds Bj 



Thence, to compute the effect (drift) on the Pioneer communication signal frequencies (uplink 
and downlink), we need only t o introduce the value of the strength of the local supercluster 
B-field: B^gc ~ 10"^ - 10"^ G |Blasi & 01into||l999| . This involves Q > 0. But then, since 



relation 51 1 implies Q < 1, and the NLED theory must be such that Lp < (hence, f + 



4ji7„ > 0) for the energy density of the EM field to be positive definite [see Novello et al.| j200"0| , 
appendix B], which entails B> B^, one can verify that the equation l |59) implies a blueshift. 
Meanwhile, the impressive accordance of the data from Voyager 1 and 2 magnetometers with 
Parker's theory (see Fig. 2) constraints Blsc to be less than 0.022 }iG within the solar system 
up to the heliopause. Hence, we may conclude that 0.01 }iG < Blsc < 0.022 within the 
solar system. By recalling that the uplink frequency of Pioneer 10 and 11 spacecrafts is v = 
2.2 GHz, one obtains for the median value Blsc = 0.018 }iG (both expressions are normalized 



70 km s-i Mpc- 



). Then Eq.i 60 ■ renders 



- I = 2.8 X 10 

V 



S-\ 



or equivalently 



dAi/ 



= 6 X 10 



-9 



Hz 



(61) 



with Av being the frequency drift pointed out earlier. 

A digression on interplanetary magnetic field and NLED effects: 

It has been pointed out that the strength of the IPMF could severly minimize the NLED effects 
because it will overrun the interstellar or intergalactic magnetic fields at heliocentric distances. 
Notwithstanding, the actual data from Voyager 1/2 spacecrafts of the IPMF average strength 
(see Fig. 2) are both consistent with a non-zero local supercluster magnetic field (LSCMF) 
amounting up to 0.022 }iG |NASA Voyager 1| jl978| ; [NASA Voyager 2 \1979,) (the accuracy 
of the measurements performed by Pioneer 10/11 magnetometers is at best 0.15 }iG, and 



^^Nonetheless, we stress that a conclusive method of fixing 7 should benefit of a dedicated laboratory 
experiment, in the same spirit that it was done, for instance, to fix the electron charge through Millikan's 
experiment 



etal 



al. 



A clean estimate of Bj from our definition of 7,, = —(B„c)^'"+^', below Eq.j59|, and the one in Novello 
{2004 : 7 = -rj^ On account that = (1 + Zc)"^ and7 = -S^ = -(Bic)*, Eq.(13) of Novello et 
|2004) rewrites B\/)Iq = lAQp^c"^ (Al). Thus fl,- = (B^cVSS^**)!''** yields Bj = ^-^'^il + ZcfBo{A2). 
Then, combining both relations (Al) and (A2) one gets : Bq = 0.02 (1 + z^)"* }iG (A3), Bi = 0.016 (1 + 
Zc)^^ }iG (A4). Now, Riess etal. found evidence for a transition from cosmic deceleration to acceleration 
at redshift Zc = 0.46 ± 0.13 [A. G. Riess etal., ApJ 607, 665 (2004)]. Inserting the latter figures in relations 
(A3) and (A4) yields : Bq = (0.005 ± 0.002) }iG (A5), Bj = (0.008 ± 0.002) ;<G (A6). Since CMB is pure 
radiation (i. e., E = Be not equal to zero on average), we consider that relations (A4) and (A6) give a better 
estimate of Bj than the one put forward in ^Novello et ar|j2004| . 
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0.022 piG for the low field system of Voyager 1/2 magnetometers Solar Wind (1980)). Besides, 
it is just beyond the Saturn orbit, ~ 10 Astronomical Units (AU), that the anomaly begins to 
be clearly observed. Surprisingly, it is just after passing the Saturn orbit that the strength 
of the magnetic field vehiculated by the solar wind gets down the strength associated to the 
insterstellar and intergalactic magnetic fields, as one can verify by perusing on Fig. 2 (see 
Refs jNASA Voyager l|^1978| ; [NASA VoyageT2| ^1978| ). Thus, since a magnetic field cannot 



shield (or block) any another magnetic field (the stronger field can only reroute the weaker 
field, otherwise it would violate Maxwell's laws), then it follows that the LSCMF has its 
magnetic influence extended upto nearly the location of the Saturn orbit, and in this way it 
forces the photons being emitted by the Pioneer spacecrafts from larger heliocentric distances 
to get accelerated due to the NLED effects. Besides, notice that |Mitra et al. l |2003 also shows 



that the local cloud of interstellar gas in HII regions does not keep out the Galactic magnetic 
field. 

"En passant", we call to the reader's attention the fact that some workers in the field have 
claimed that the effect should have shown up already at the small distance corresponding 
to Mars, Jupiter or Saturn orbits, because of the high technology involved in the tracking 
of planet orbiting spacecrafts as Galileo and Cassini or the Mars' nonroving landers, which 
would allow to single out the anomaly at those heliocentric distances. However, as those 
spacecrafts are inside the region where the solar wind dominates, this definitely precludes 
the NLED photon acceleration effect to show up at those distances since the much higher 
magnetic field there would introduce a negligible NLED effect, and as stated below the solar 
pressure influence on the signal frequency is still large. 

On the other hand, although one can use the time of flight of photons during tracking of 
planets with orbiting spacecrafts (by combining range and Doppler data over a spacecraft 
orbit) to tightly determine the range from Earth to that given planet's center of mass, the 
impediment to single out the radio-signal frequency shift remains the same pictured above: 
from one side the strength of both the host planet magnetic field and the solar magnetic field 
at those distances are still large, what blocks the nonlinear action of the LSCMF, and from the 
large up to 20 AU so as to allow the show up of the NLED frequency shift which is much 
smaller. Moreover, within a heliocentric distance ~ 100 AU the IPMF keeps stronger. Thus, 
it reduces for all practical purposes the IPMF contribution to the effects of NLED (see further 
arguments from our direct estimate of Bj in Footnote 13 above, and in Ri ess etaT] j2004) ). 



leaving room for the sole contribution of the residual IGMF in the solar system. 
Finally, the new frequency shift that is predicted by NLED is not seen yet in the laboratory 
because of some of the following reasons: a) the most important, the strength of the Earth 
magnetic field is much larger than the one required in the NLED explanation of the anomaly 
for the effect to show up, and b) the coherence time t = 1/Av of EM waves in present 
atomic clocks (frequency width Av > 0.01 Hz, or otherwise stated ct < 0.2 AU) is too short 
as compared to the time of flight of photons from Pioneer 10/11 spacecrafts past 20 AU. 
Nonetheless, if the conditions demanded by our model were satisfied this effect will certainly 
be disentangled in a dedicated experiment where, for instance, the Earth magnetic field is 
kept outside the case containing an experimental set up where a very weak magnetic field is 
maintained inside, a source of photons set to travel and a receiver data-collecting. 
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Fig. 2. Data from Voyager 1/2 spacecrafts of the interplanetary magnetic field (IPMF) 
average strength. Subfigure 1 presents measurements by Voyager 1 of the strength of the 
IPMF [nT] as a function of time [Yr]. The continuous line represents the predictions of 
Parker's model, and the dots the Voyager 1 data. Subfigure 2 shows the dependence with the 
distance [AU] of the IPMF [nT], as detected by Voyager 2, against the theoretical prediction. 
Credit: http:/ /spacephysics.ucr.edu/images/swq2-04.jpg, and 
http: / / ir\terstelIar.jpI.nasa.gov/interstellar/ probe / interstellar /images/ 
07BupstrmSuess-25Ig.gif 
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and for extragal actic clusters as Abell 400, Hydra A measurements provide B^g-c ' 



3 - 20 X 10-^ G 
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Following Hadamard |Hadamard| jl903| , the surface of discontinuity of the EM field is 
denoted by E. The field is continuous when crossing E, while its first derivative 
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7. APPENDIX: The method of effective geometry 

Following Hadamard (1903), the surface of discontinuit)j^ of the EM field is denoted by 
E. The field is continuous when crossing E, while its first derivative presents a finite 
discontinuity. These properties are specified as follows 




where the symbol 

[ffiy]j; = lim^-^o+(/l2:+(5 - J\i.-s) 
represents the discontinuity of the arbitrary function / through the surface E 
is called the discontinuity of the field, and k,\ = 3,\E is the propagation vector 
symbol "|" stands for partial derivative; 

^*Of course, the entire discussion onwards could alternatively be reptirased in terms of concepts 
more familiar to the astronomy community as that of light rays used for describing the propagation of 
electromagnetic waves in the geometric optics approximation. 



(63) 

. The tensor f^y 
. In Eq. |62), the 
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Here-after we want to investigate the effects of nonlinearities of very strong magnetic fields 
in the evolution of electromagnetic waves described onwards as the surface of discontinuity 
of the electromagnetic field (represented here-to-fore by Ffn,). For this reason we will restrict 
our analisys to the simple class of gauge invariant Lagrangians defined by L = L{F). From the 
least action principle we obtain the following field equation 

(Lfp-)|l^, = 0. (64) 

Applying the Hadamard conditions (j62| and (j63f to the discontinuity of the field in Eq.|64} we 
obtain 

Li:p"'h,+2LFF^Ff"-'ky = 0, (65) 

where f is defined by f = f^f^g. Both, the discontinuity conditions and the electromagnetic 
field tensor cyclic identity lead to the following dynamical relation 

flivk^ + fyAk,, + hjiky = . (66) 
In the particular case of a polarization such that f = 0, it follows from Eq.j64) that /''"fc,, = 0. 



Thus using this result, and multiplying Eq.i 66 by k'^ we obtain 



UvkV = . (67) 

This equation states that for this particular polarization the discontinuity propagates with 
the metric fpiv of the backgroimd space-time. For the general case, when f 7^ 0, we multiply 
Eq.l|66| by k^^^Ft'" to obtain 

?/c,fc,,^f'^' + 2F''7,^fcAfc,,=0. (68) 

From this relation and Eq.|65) we obtain the propagation law for the field discontinuities, in 
this case given as 

Lf^''^-4LFfFl'f'''')fc,,fcv = 0, (69) 



where F^F"^ = —B h}''" — B''B''. Eq.i 69 ■ allows to interpret the term inside the parenthesis 



multiplying fcf'fc^ as an effective geometry 

^eff = ^f^'"-4i^fFf'f"- (70) 
Hence, one concludes that the discontinuities will follow geodesies in this effective metric. 

& APPENDIX-1: NPS theory applied to cosmology 

To discuss the evolution of a universe model driven by the NPS NEED, the electromagnetic 
( EM) field described by Eq.pi} can be taken as source in Einstein equations to obtain a toy 
model for the evolution of the universe which displays accelerate expansion. Such phase of 
acceleration runs into action when the nonlinear EM term takes over the term describing other 
matter fields. This NEED theory yields ordinary radiation plus a dark energy component 
with zv < —1 (phantom-like dynamics). Introducing the notatioriM the EM field can act as 



^^Due to the isotropy of the spatial sections of the Friedman-Robertson-Walker (FRW) model, an 
average procedure is needed if electromagnetic fields are to act as a source of gravity [Tolman| (l934[ . 
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a source for the FRW model if (£,)|^, = 0, (B,)|^, = 0, (£,■B;)|^, = 0, (£,Ey)|^ = - ^E^gjj, and 
{BiBj)^^ = "jB^^/y]^ When these conditions are fulfilled, a general nonlinear Lagrangian 
L(F) yields the energy-momentum tensor (Lf = dL/dF, Lff = ij^L/dF^|^ 



{Tf>v)\^ = {p+p)v^,Vy-pg^,y, (71) 

p = -L-4E^Lf, p = L+ ^(£2-2B2)Lf, 



Hence, when there is only a magnetic field, the fluid can be thought of as composed of a 
collection of non-interacting fluids indexed by k, each of which obeys the equation of state 

Pj. = (t ~ -'^) Pt<' composed of ordinary radiation with pj = j and of another fluid with 

equation of state P2 = — 3 P2- It is precisely this component with negative pressure that 



may drive accelerate expansion through Friedmarm equations, as was shown in Novello et 
|al]|2004| 



Thus a volumetric spatial average of a quantity X at the time f by = limi/^vi) y [ X^— g d^x, 

where V = J ^—g (Px , and Vq is a sufficiently large time-dependent three-volume. (Here the metric 
sign convention (H ) applies). 

^^Let us remark that since we are assuming that (B,) =0, the background magnetic fields induce no 
directional effects in the sky, in accordance with the symmetries of the standard cosmological model. 

^''Under the same assumptions, the EM field associate to Maxwell Lagrangian generates the 
stress-energy tensor defined by Eq.j7l| but now p 



